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We propose a new linearizable model for the nonlinear photocurrent-voltage characteristics of 
nanocrystalline Ti02 dye sensitized solar cells based on first principles and report predicted values 
for fill factors. Upon renormalization diverse experimental photocurrent-voltage data collapse onto 
a single universal function. These advances allow the estimation of the complete current-voltage 
curve and the fill factor from any three experimental data points, e.g., the open circuit voltage, 
the short circuit current and one intermediate measurement. The theoretical underpinning provides 
insight into the physical mechanisms responsible for the remarkably large fill factors as well as their 
known dependence on the open circuit voltage. 

PACS numbers: 84.60.Jt, 72.80.Le, 05.40.-a, 85.60.-q 



I. INTRODUCTION 

An important feature of photovoltaic solar cells and 
of diverse optoelectronic devices studied in semiconduc- 
tor physics concerns their current-voltage characteris- 
tics [3, 0, 0, 0|- The pioneering work that led to the 
invention of Gratzel or dye sensitized solar cells became 
a milestone in the study of photovoltaic devices @, @. 
Previous theoretical and experimental studies have iden- 
tified the dependence of the photocurrent and photovolt- 
age on radiant power [7| , but not the precise nonlinear de- 
pendence of the photocurrent on the photovoltage under 
conditions of constant radiant power. Moreover, vari- 
ability in the manufacturing process of dye sensitized 
solar cells can lead to differences — e.g., variables in- 
clude the choice of dye, the sintering temperature, thick- 
ness of the nanoporous Ti02 film and choice of chemical 
treatments. This diversity leads to significant qualita- 
tive and quantitative variation in photocurrent-voltage 
characteristics and of the relevant quantities such as the 
open circuit voltage V oc or the fill factor. Such variability 
has discouraged attempts to identify (possibly "hidden" ) 
dynamical patterns that could yield important insights 
into the regenerative photoelectrochemical mechanisms 
that underlie the conversion process. Given the variabil- 
ity and diversity in the characteristics, which properties 
remain universal and which nonuniversal? More impor- 
tantly from an experimental point of view, how can we 
quantitatively model the photocurrent-voltage character- 
istics, based on fundamental theoretical principles? Here 
we answer these questions by deriving from first princi- 
ples an analytical expression for the photocurrent. 



The topic of solar energy in general |8 



and 

dy e sensitized solar cells in particular [a Ifil; LZI 
EE EE EE EH attracts broad interest from diverse sectors 
of society, due to technological, economic, political and 
environmental considerations. The growing scientific in- 
terest in dye sensitized Ti02 solar cells stems from their 
unusual features and mode of operation that distinguish 
them from Si solar cells: (i) efficient charge separation 



due to ultra-fast injection of electron from the dye on 
picosecond and subpicosecond time scales [H], E3> ( u ) 
conduction consisting only of injected electrons rather 
than electron- hole pairs [5|, [6| , due to the wide bandgap 
of the semiconductor T1O2; (ii) high optical density due 
to the extremely large surface area of the dye sensitized 
nanoporous semiconductor [l4| ; (iv) negligible charge re- 
combination with the oxidized sensitizer dye 
and (v) very high quantum yields One important 

fact that will contribute towards derivation and subse- 
quent interpretation of the current-voltage characteris- 
tics concerns how the experimentally measured recombi- 
nation current density vanishes at short circuit Q — in- 
dicating that the only significant recombination pathway 
proceeds via back electron transfer into the electrolyte. 
Indeed, charge recombination between redox species (IJ 
ions) in the electrolyte and conduction band electrons 
localized at the nanoporous interface result in subopti- 
mal photovoltage levels — thus limiting the conversion 
efficiency [T^, Ho, - 



II. METHODS 

We begin by assuming that the number of electrons 
injected into the conduction band depends only on the 
incident radiant power — in fact the known very high 
quantum yields justify this assumption. This assumption 
allows us to express the recombination current density J r 
as a function of the photocurrent density J and the in- 
jection current density. Since J r vanishes as short circuit, 
the injection current equals the short circuit current J sc , 
so that 



J 



(1) 



For a well mixed solution with identical surface and bulk 
concentrations (typical for small current densities), the 
Butler- Volmer equation [13, [H| leads to the following 
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Figure 1: Typical photocurrent voltage characteristics under 
a radiant power of 1.5 AM for dye sensitized solar cells, taken 
from ref. The four data sets have different characteris- 
tics due to varied chemical treatments (see text), yet in each 
case the theoretical curves (solid lines) corresponding to Eq.|6] 
can account well for the experimental curves. Inset shows ap- 
proximate logarithmic relation for the recombination current 
density J r versus voltage, not inconsistent with Eq. [3] 



expression for the recombination current density J r 



Jr = Jo 



exp(-a c M/r?) - exp(a A ufrj) 



(2) 



where Jo denotes the exchange current density, u the 
number of electrons transferred in the reaction (and con- 
sequently, the order of the rate of reaction for recom- 
bination for electrons) a A and a c the anodic and ca- 
thodic transfer coefficients, T) = V the overpotential and 
/ = q/k B T. With some simplification [tJ, this equation 
becomes 



J r = qk ot c m n% a [exp(uaqV/k B T) - 1] 



(3) 



where a = a A and fc Gt denotes the back electron transfer 
rate constant, c the concentration and m the order of the 
reaction for the oxidized species, q the electronic charge 
and no represents value in dark conditions of the electron 
population in the semiconductor. 

How do the prefactors depend on voltage V? The back 
electron transfer rate constant varies with radiant power 
and the maximum, i.e. open circuit, photovoltage V oc , 
however we do not expect it to depend on the photovolt- 
age for fixed radiant power. The Nernst Equation for 
the potential in terms of the concentrations of the oxi- 
dized and reduced species holds valid only under equi- 
librium conditions, yet we know that c varies across the 
electrolyte. Indeed, since the reduced species greatly ex- 
ceeds the oxidized species, we can safely conclude that 
the voltage varies as AV(x) w — (k B T/q) \n(c°/c(x)), 
where c(x) denotes the concentration at a position x 
across the cell (i.e., electrolyte), c° denotes a reference 
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Figure 2: Renormalized photocurrent J^ R ' versus V/V c for 
the curves shown in Fig. [TJ The data collapse onto a single 
universal function stable over a variety of chemical treatments 
and solar cell variability. We have chosen v s = 1/40, based on 
the value of k B T at room temperature, for illustration, cor- 
responding to our estimate of an upper bound for fill factors, 
however we could have renormalized the photocurrent to any 
fill factor. Indeed, we can linearize the curves, as shown in 
the inset and explained in the text. The dashed line traces 
the upper and lower uncertainties corresponding to an error 
of 1% in the short circuit current. Notice how remarkably the 
data collapse onto a straight line, within the error tolerance. 



(or mean) concentration. Nonetheless, we still do not 
know, a priori, exactly how it varies with the potential 
at the semiconductor-electrolyte interface as the external 
load (i.e., impedance) varies, because of the out of equi- 
librium conditions. In this context, one important clue 
comes from the dependence on the photovoltage on the 
electron population. Electrons act as charge carriers in 
the T1O2 — just as the redox species do in the electrolyte. 
In the semiconductor, injected electrons shift the Fermi 
level, so that n = hq exp(qV/k B T). This exponential 
dependence on potential, together with the exponential 
dependence of c on voltage across the electrolyte, hints 
at a similar, i.e. exponential, dependence of c on V at 
the interface as the external load varies: 



V = V oc - (jk B T/q) ln(c oc /c) 



(4) 



Here c oc represents the concentration under open cir- 
cuit conditions and 7 represents a free parameter in the 
model, quantifying the fraction of the voltage variation 
that affects the oxidized species concentration. Since the 
triiodide concentration cannot vary very much in the liq- 
uid, we cannot expect small 7 close to equal unity since 
this would imply can. For now we only mention that 
one would naively expect a positive 7 since electron in- 
jection and dye regeneration (associated with larger V) 
produce oxidized species. 
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Figure 3: (a) Theoretically predicted fill factors FF and ratio 
Vopt/K) C of optimal operating voltage Vopt to open circuit 
voltage V oc versus V oc /Vs, where Vs represents a characteristic 
voltage (Eq.Q. We estimate a worst-case lower bound on the 
fill factor using Voc/V s = 1, and an idealized upper bound 
using Voc/Vs = 40. Inset shows a more complete picture, (b) 
Comparison of theoretical and experimental values for FF for 
the data shown in Fig[T] Useful real cells will likely have FF 
within the range shown. 



III. RESULTS 

These considerations immediately lead to an analytical 
expression for the photocurrent J as a function of the 
voltage V across the cell: 



1 



exp(mq(V - V oc )/-/k B T) 
exp(uaqV oc /k B T) - 1 

exp(uaqV / k B T) — 1 



(5) 



Here m « 2 because of the second order reaction and 
ua « 0.7 @, EH- This analytical expression must hold 
true for all dye sensitized solar cells, yet in the context of 
providing greater clarity and insight we can render it sim- 
pler but more useful. Moreover, from a practical point of 
view, we can simplify it further by making additional yet 
realistic assumptions. Eq. [3] for the recombination has 



validity in the large voltage (V > 80mV) regime fA [l3l|. 
On the other hand, below this potential, the recombina- 
tion current becomes negligible and uninteresting. For 
any useful cell, n 3> uq by many orders of magnitude, so 
that we can reasonably well approximate Eq. [5] with 

~exp(V oc /V s ) - exp(V/V s 



J J a. 



where the potential 
V s = {k B T/q)- 



exp(V oc /V s ) - 1 
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-Volts 



(G) 



(7) 



ua + m/j 40(0.7 + 2/7) 

represents a characteristic scale of the exponential decay. 
Specifically, V s quantifies the photovoltage drop corre- 
sponding to a decrease in recombination current density 
by a factor of 1/e where e here denotes Euler's number. 
Eqs. [5] and O represent the first out of three new results 
of this article. 

Fig. [T] compares the model with photocurrent-voltage 
curves taken from ref. Q, of untreated and pyridine 
derivative-treated [RuL2(NCL)2]-coated nanocrystalline 
Ti0 2 electrodes in CH 3 CN/MNO (50:50 wt %) con- 
taining Li(0.3M) and l2(30mM), for a radiant power of 
100mW/cm 2 (AM 1.5). The electrodes had treatment 
with following substances: 3-vinylpyridine (VP), A-tert- 
butylpyridine (TBP), and poly(2-vinylpyridine) (PVP). 
The good agreement with the data validates the model 
represented by Eqs. [5] and [SI 

The largest possible power output divided by J sc V oc 
defines the fill factor FF. Notice that V s changes the fill 
factor (via 7). A value V s — > 00 (corresponding to purely 
resistive or Ohmic behavior) leads to FF=l/4, whereas 
V s — > leads to unity fill factor — perfect but theoreti- 
cally impossible except at T = K. Most dye sensitized 
solar cells have FF=0.6-0.7 (Fig. [J). 

We now turn our attention to the question of whether 
a single universal current-voltage relation can describe 
all TiC>2 solar cells. According to the theory presented 
above, all dye sensitized solar cells must satisfy Eq. [5] if 
not Eq. [6] If we renormalize the photovoltage to obtain 
an adimensional measure V* = V/V oc , then every single 
dye sensitized solar cell must satisfy the following relation 
for an idealized renormalizcd photocurrent: 

Vs/Voc^s 



1 - 



J<*> = 



1 — (J/ J sc ) I 1 - exp(-V r oc /y s ) 



1 - exp(-l/u s ) 



(8) 

Here v s fixes the shape or fill factor of the renormalizcd 
photocurrent. Fig. [2] shows the predicted data collapse. 
We have chosen a value v s = 1/40, due to its significance 
for an idealized solar cell with maximum FF (see below) 
at room temperature. However, we can obtain data col- 
lapse for any v s (not shown). This is perhaps more clear 
if we linearize the curves. We define coordinates 



X* = l-V/V oc 
Y* = -V s /V oc ln 



1--H 1 - cxp(-V oc /V B ) 



(9) 
(10) 
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The inset of Fig. [2] shows how the data collapse onto a 
straight line. All dye sensitized solar cells thus follow the 
same universal pattern of photocurrent- voltage behavior. 

We next consider the problem from the point of view 
of scale invariance symmetry. The fill factor cannot 
depend on J sc , since it cancels in the power ratio. 
It also remains invariant under a scale transformation 
V oc — ► ^V oc , V s —> XV S . In fact, no dilation can alter a 
ratio of geometric areas. The invariance of FF for arbi- 
trary A implies that FF can depend on V oc and V s only 
via their ratio: 

FF = FF{V oc /V s ) . (11) 

The exact functional dependence appears to involve a 
transcendental equation. We are still attempting an an- 
alytical solution using the Lambert W function. Never- 
theless, it is susceptible to numerical solution. Fig. [3]^a) 
shows FF as a function of V oc /V s . 

We next comment on the values typically found for V s 
and their physical significance. The values found cor- 
respond to negative 7 and thus suggest that the con- 
centration of redox species (LJ ions) decreases with the 
photovoltage. This may at first seem counter-intuitive. 
Indeed, higher voltage suggests larger electron popula- 
tion and more injection. Moreover, the regeneration of 
the dye creates Fg~ species, in the proportion of one ion 
for every two electrons injected. 

So do we face an apparent inconsistency? The impor- 
tant fact, mentioned earlier, of zero recombination cur- 
rent density J r = under short circuit conditions, hints 
at the correct explanation: the rate of regeneration of 
the oxidized dye depends not on the photovoltage (zero 
under short circuit) but rather on the rate of electron in- 
jection — thus only on the open circuit photovoltage, or 
alternatively, on the radiant power. The finding agrees 
with the expectation of a smaller depletion layer for more 
external current drain. 

The above findings allow us to estimate lower and up- 
per limits for FF. Purely Ohmic behavior corresponds 
to V s — * 00 and FF=l/4, however we cannot imagine 
this scenario. For any useful device, the largest conceiv- 
able value of V s should not exceed V oc , which gives us 
a lower bound for FF of FF=0.31 and an optimal op- 
erational voltage of V OPT = 0.55 V oc . By considering 
V s = 1/40 Volts, i.e. idealizing ua = 1 , m/7 = 0, and 
V oc = 1, we arrive at an upper bound of FF=0.88 and 
Vopt = 0.91 Voc, as shown in Fig. Ela). For ua = 0.7 
we obtain slightly smaller FF. Fig. [3{b) allows one to es- 
timate one among V s , V oc and FF from the other two 



and will thus find practical application. We estimated 
the error bars for the experimental points from the re- 
gression fits used to arrive at the values of V s . Devices 
that we constructed locally had values of FF within these 
bounds. 

Finally, our findings explain the very large fill factors 
of dye solar cells. The recombination current becomes 
insignificant as soon as the voltage drops to V = V oc — V s 
(Eqs. [3J El). If V s <C V oc (as in fact happens), then the 
photocurrent jumps from zero to close to its short circuit 
value even if the voltage only drops slightly (i.e., by V s ). 
Notice from Fig. [3] that increases in V oc — e.g., due to 
greater radiant power — should indeed lead to higher FF 
if V s varies much less than V oc , as in fact occurs. 



IV. DISCUSSION AND CONCLUSION 

We briefly comment on conversion efficiency. Since we 
cannot expect to change V s significantly, further increases 
in efficiency will depend mainly on increasing the open 
circuit photovoltage, which in turn also requires reducing 
losses due to charge recombination at the nanoporous 
interface — the main challenge indeed. Moreover, we 
know that c ~ n 1 / 7 , where 7 has the negative value 7 w 
—6. In other words, we find evidence of a fractional power 
law or scaling exponent, indicating self-affine behavior. 
We hypothesize that the negative value arises due to the 
fact that larger n leads to greater recombination, which 
consumes the oxidized species. Localization effects and 
and transport properties play an important role in this 
context [181 ]. 

In summary, our theoretical results appear to account 
well for the observed behavior of real dye sensitized solar 
cells and seem to provide new insights into their function- 
ing. Among the important results reported here, we note 
that Fig. [3] allows one to calculate FF knowing V oc /V s or 
vice versa. Moreover, knowing either one or the other, 
one can readily obtain the entire photocurrent-voltage 
curve, via Eq. [6j From just 3 points in the photocurrent- 
voltage curve, one can reconstruct the entire curve. The 
findings reported here may thus allow further advances 
and eventually lead to technological innovations. 
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